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ABSTRACT 

We propose to determine the mass and the radius of a neutron star (NS) using three measurable 
mass-radius relationships, namely the "apparent" radius inferred from neutron star thermal 
emission, the gravitational redshift inferred from the absorption lines, as well as the averaged 
stellar mass density inferred from the orbital Keplerian frequency derived from the kilohertz 
quasi periodic oscillation (kHz QPO) data. We apply the method to constrain the NS mass and 
the radius of the X-ray sources, IE 1207.4-5209, Aql X-1 and EXO 0748-676. 
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1 INTRODUCTION 



Measuring the mass (A/) and the radius 7? of a neutron star (NS) 
is very important in both nuclear physics and gravitational physics 
since it allows us to constrain the NS matter compositions (e.g. 
neutrons or quarks) at super-density (see, e.g. Lattimer & Prakash 
2004; Miller 2002; Li et al. 1999; Cheng et al. 1998; Haensel et al. 
1986; Alcock et al. 1986), and to trace the particle motion behav- 
ior in super-strong gravitational field (see, e.g. van der Klis 2000, 
2006). 

NS masses have been measured in the radio pulsar binary sys- 
tems with high precision (see, e.g. Kaspi et al. 1994; Bailes et al. 
2003). For example, the masses of the two NSs of the Hulse-Taylor 
pulsar system, PSR 1913+16, have been measured as M=1.41Mq 
and M=1.38 Mq, respectively (see, e.g., Manchester & Taylor 
1977); and those of the recently discovered double pulsar sys- 
tem, PSR J0737-3039, are M=1.34M0 and M=1.25 Mq, respec- 
tively (see, e.g., Lyne et al. 2004). The mass constraints of 26 
neutron stars in the binary radio pulsar systems have been derived 
by Thorsett & Chakrabarty (1999), who present a remarkably nar- 
row underlying gaussian mass distribution at M=1.35±0.04Mq . 
In the low mass X-ray binaries (LMXBs) a lower mass limit 
M > 0.97 ± 0.24Mq was derived for the NS in 2A 1822-371 
(Jonker et al 2003), and the NS mass in Cyg X-2 was derived as 
M = 1.78 ± O.23M0 (Orosz & Kuulkers 1999). 

In contrast, there is no method so far to directly measure the 
radius of a NS. A conventional value of NS radius is about 15 km, 
which is derived by requiring that the averaged dipolar magnetic 
field strength (i.e. ~ 10^^ G) inferred from radio pulsar spindown 



(see, e.g. Manchester & Taylor 1977) is consistent with the one 
implied by the cyclotron line emission in some X-ray binaries (see, 
e.g. Makishima et al. 1999; Truemper et al. 1978). The uncertainties 
in determining NS radius make it difficult to accurately measure 
the NS equation-of-state (EOS). The detail of matter composition 
inside a NS has been an open issue. 

With the launches of X-ray space observatories, especially the 
recent XMM-Newton, Chandra, and RXTE (see, e.g. Jansen et al. 
2001; van der Klis 2000, 2006), some M - R relationships have 
been appropriately measured or estimated. These include (1) the 
"apparent radius" Rao estimated from the thermal emission of iso- 
lated NSs (e.g. Geminga, Caraveo et al. 2004; IE 1207.4-5209, 
Bignami et al. 2004; RX J 1856.5-3754, Truemper et al. 2004; Bur- 
witz et al. 2003; Burwitz et al. 2001) or NSs in binary systems 
(e.g. Aql X-1, Rudledge et al. 2001); (2) the gravitational redshift 
implied from the absorbtion lines in EXO 0748-676 and in 4U 
1700-F24 (e.g., Cottam et al. 2002; Tiengo et al 2005) and (3) the 
averaged NS mass density inferred from the orbital Keplerian fre- 
quency used to interpret the kHz QPOs discovered by RXTE (see, 
e.g., Miller et al. 1998; van der Klis 2000, 2006; Zhang 2004). 

Generally speaking, if two of the three above M — R rela- 
tions are well "measured", one can solve both M and R indepen- 
dently. This is the main topic of this paper. We will discuss the 
methods in detail and apply them to several sources. Once NS M 
and R known, we can evaluate the matter compositions of star 
by comparing with the representative EOSs, as shown in the 
three figures we chosen the EOSs of strange matter (CSl and 
CS2), normal neutron matter (CNl and CN2) and pion conden- 
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sate in the star core (CPC) (see, e.g. Lattimer & Prakash 2001; 
Cook et al. 1994). 



MEASURING NS MASS AND RADIUS WITH M 
RELATIONS 



R 



(1) The "apparent radius" Rao inferred from the thermal emis- 
sion data is only an upper limit of the NS radius. It is related to 
the "true" radius R through (see, e.g., Haensel 2001; Lattimer & 
Prakash 2001; Thome 1977), 



Rao = R/^/l-Rs/R 



(1) 



where Rs — 2GM is the Schwarzschild radius. From Eq.O, 
the maximum value of the mass can be easily obtained as 
M = 1.3(i?oo/lO^cm)M0 when R = 0.58i?oo. We stress that 
the apparent radius is not a directly observable quantity al- 
though derivable (in principle) from two observable quantities, 
namely the "apparent blackbody luminosity" (the blackbody 
luminosity seen by a distant observer) and the "apparent black- 
body temperature". Both quantities are derived by fitting X-ray 
data under the assumption of perfect blackbody, which is, in 
general, not true (see, Haensel 2001). In particular the inclu- 
sion of an atmosphere above the NS crust can vary the estimate 
of the NS apparent radius by more than a factor 2 (see Haensel 
2001, section 3.1), thus the 'measured' apparent radius is still 
model dependent. 

(2) The gravitational redshift near the neutron star surface pro- 
vides another M-R relation, giving the ratio between the NS mass 
and the NS radius (e.g. Miller 2002). The approach becomes rel- 
evant with the recent discovery of the absorption features in EXO 
0748-676 (Cottam et al. 2002). According to the definition of the 
gravitation redshift 2 = A A/ A ~ GAf/7? (for a weak gravitational 
field, where A is the wavelength of the emission), for a non-rotating 
spherical gravitational source in the Schwarzschild spacetime, we 
get (see, e.g., Sanwal et al. 2002), 



1 + z = 1/^/1 ~Rs/R., 



or, in terms of M/R, 



m/Re, = f{z) = 



z{l + z/2) 

0.15(1+2)2 



(2) 



(3) 



where Re = R/lO^cm, and m — M/M© is the mass in unit of 
solar mass. In the case of weak gravitational fields (z <^ 1), f{z) ~ 
2;/0.15 is approximately obtained. 

(3) A third M — R relation can be obtained by analyzing the 
kilohertz quasi periodic oscillation (kHz QPO) data. This field has 
been greatly benefited from the RXTE mission, which discovered 
the twin (upper and lower) kHz QPOs or millisecond variability in 
the X-ray data of about 20 NS LMXB systems. The upper kHz QPO 
frequency is generally interpreted as the Keplerian frequency uk of 
orbital materials at some preferred radius r near the NS surface (van 
der Klis 2000, 2006; Morsink 2000; Stella & Vietri 1999; Zhang 
2004), i.e.. 



UK = 



GM 



= 1850 (Hz)AX^''^ , X = R/r 



(4) 



with the parameter A = (m/i?!)^^^ 

or Re, = 1.27m^''^(A/0. 7)^2/3^ ^jjgj.g ^j, jj^g j^g ^.^(Jj^g -^^ 

unit of 10 km. The physical meaning of the parameter A is such 



that the quantity represents the "measurement" of the averaged 
NS mass density, namely. 



M/i?^ ~ 10''*(g/cm^) (4/0.7)^ 



(5) 



Meanwhile, by interpreting the twin kHz QPOs as the Keple- 
rian motion and the Alfven wave oscillation mechanism (see, e.g., 
Zhang 2004), the averaged mass density parameter A can be in- 
ferred from the simultaneously detected twin kHz QPO frequen- 
cies. For typical kHz QPO sources, yl ~ 0.7 is approximately ob- 
tained. If, instead, only one single kHz QPO frequency is detected, 
one can not derive the exact value of A from the model, but can esti- 
mate the lower limit of A and the upper limit of M according to the 
following arguments, as proposed by Miller et al. (1998). In order 
to explain the observed X-ray flux saturation at the twin kHz QPO 
frequencies (e.g. van der Klis 2000, 2006; Zhang et al. 1998), it is 
usually believed that the maximum Keplerian frequency vk occurs 
either at the star surface {r = R: X = 1) or at the innermost stable 
circular orbit (ISCO) with the radius -Risco = 3i?s = 6GM. The 
two inequalities are therefore derived from Eq.lO of Miller et al 
(1998), i.e. 



and 



A > //jf/1850(Hz) 



m < 2200(Hz)/i^K 



(6) 



(7) 



In addition, Burderi & King (1998) have derived & M — R relation 
for the millisecond accretion-powered X-ray pulsar SAX J 1808.4- 
3658 with the measured spin frequency 401 Hz through comparing 
the stellar radius with its co-rotation radius of the magnetosphere- 
disk. They inferred a value of A> 0.62. 

To summarize, hitherto there are roughly three kinds of well- 
known M — R relationships that could be measured or estimated 
from the data, as described above. Knowing two of the three in a 
system would in principle lead to the determination of both M and 
R. In reality, errors are large for these relationships. Knowing all 
three relations for a source would greatly diminish the errors. As 
an illustration, below we will use the method to constrain M and 
R for several sources. 



2.1 Constraining M and R with Rao and z: the case of IE 
1207.4-5209 

If the "apparent radius" Rao and the gravitational redshift z are 
known, both M and R can be calculated by the following two for- 
mulae, derived from Eqs.Q and (O, 



Ra = Rao(>/(l + z) , 

m = f{z)Rao(i/{l + z) 



(8) 



(9) 



withiiooe = Rao/lO^cm. 

In the spectral analysis of IE 1207.4-5209, Bignami et al. 
(2004) found that the best fitting continuum model includes two 
blackbody components, i.e. a cooler component with a temperature 
KT = 0.163±0.003keV and an emitting radius of i?e = 4.6±0.1 
km, and a hotter component with KT = 0.319 ± 0.002 keV and 
Re = 0.83±0.03 km. If we take the cooler component as the emis- 
sion from the whole star, then the cooler emitting radius 4.6±0. 1 
km can be regarded as a measured lower limit of the "apparent ra- 
dius" Rao- In addition, an emitting radius R=4.6 ± 0.1(d/2fcpc) 
km in IE 1207.4-5209 is obtained by the best fitting continuum 
model of the blackbody function and by the distance of 2kpc 
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(Bignami et al. 2004), however as claimed by Haensel (2001) 
that an emitting radius may be doubled if the atmosphere 
model of the photon spectrum is considered, i.e. the emitting 
radius may be as high as ~9.4 km. 

On the other hand, two absorption features in the spectrum of 
IE 1207.4-5209 were detected with the Chandra ACIS (Sanwal et 
al. 2002). Interpreting these lines as Helium atomic lines, a grav- 
itational redshift is inferred as 2; ~ (0.12 — 0.23) (Sanwal et al. 
2002). This conclusion is, however, inconclusive, since the lines 
may be due to the cyclotron origin (e.g. De Luca et al. 2004). Here 
we take this gravitational redshift measurement as an unconfirmed, 
tentative case. 

In Figure [T] we present the M — R diagram of NSs with 
various EOSs and present the above two M — R constraints in 
the same plot. For the possible parameter ranges of R^a and 
z = 0.12 - 0.23, i.e. Roa = 4.5 - 4.7 (km) as estimated from 
the blackbody assumption (Bignami et al 2004), Roc = 9.4 (km) 
while the atmosphere model of photon spectrum is taken into 
account (Haensel 2001), and Raa ~ 20 (km) assumed under 
the condition that the uncertainty in the distance of the source 
1E1207.4-5209 to the earth is considered (Truemper 2005), we 
constrain the parameter space within the area enclosed by the 
lines labeled Rl, R2, R3, and those labeled z=0.12, 0.23, where 
the implied M —7? ranges are given. In Figure[ll we also plot the 
lines ot R = Rs = 2GM (the solid straight line labeled as Rs, 
below which the M — R values will be permitted for NSs) and 
R ~ i?isco ~ 6GM (the middle straight dotted line labeled as 
ISCO). The latter corresponds z — 0.23 (this value is fortu- 
itous coincidence with the observationally inferred redshift up- 
per limit of 1E1207.4-5209, see Sanwal et al. 2002) and divides 
the M — R diagram into the R < iiisco region (above the 
line) and the R > i?isco region (below the line). From Fig- 
ure [1] the inferred M - R value regime for 1E1207.4-5209 fa- 
vors the existence of the exotic matter inside NS. However, it 
is pointed out by Truemper (2005) that the distance d from IE 
1207.4-5209 to the earth is not so certain and d = 1.3 ~ 3.9fcpc 
has been given (see, e.g. Pavlov et al. 2004). If we consider 
the atmosphere model of photon spectrum and the distance 
d = 3.9kpc, then the derived "apparent radius" can be as large 
as ISkm (4.6*2*3.9/2=18), which will weaken the conclusion of 
IE 1207.4-5209 including a strange star because the inferred 
M ~ R range in Figure[T]is enlarged to cover the EOSs of nor- 
mal neutrons. 

If one assumes the validity of the above measured R^c (4.5 
- 4.7 km) and z (0.12 - 0.23), then one gets M = 0.34 ± O.OQMq 
and ii = 4.2 ± 0.1 km for IE 1207.4-5209 according to Eqs.© 
and @. In this case, one can conclude the NS EOS of IE 
1207.4-5209 infers the composition of strange quark matters 
(e.g. Xu 2005). Apparently, this conclusion is model dependent 
and based on the blackbody assumption and the distance of 
2kpc. 



2.2 Constraining M and 7? with Roc and A: the case of Aql 
X-1 

If the apparent radius Roc and the parameter A are known, both 
Al and R can be calculated by the following two formulae, derived 
from Eqs.O and 




Re = i?oo6/\/l + 0.15(71/0.7)2/?^ 



(10) 
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Figure 1. The M — R diagram. Five representative EOSs are shown 
(see also, e.g. Miller 2002): stars containing strange quark matter (CSl 
and CS2), stars made of normal neutron matter (CNl and CN2), and 
stars with pion condensate cores (CPC) (Lattimer & Prakash 2001; 
Cook et al. 1994). The straight lines are constant z lines: the solid 
(Rs) and the dotted lines (ISCO) represent that the stellar radius is 
one and three times of the Schwarzschild radius, i.e. R = 2GAI (be- 
low the line Rs the NS M — R. regimes are possible, see e.g. Haensel 
et al. 1986) and R = 6GM, respectively. The constant redshift lines 
z = 0.12 and z = 0.23 are also plotted. The concave-down parabo- 
las are constant Roc lines: Roc = 4.5, 9.4, 18 km labeled by Rl, R2 
and R3 correspond the lower limit of the apparent radius of 1E1207.4- 
5290 (4.5=4.6-0.1) inferred by the blackbody assumption (Bignami et 
al. 2004), (9.4=2*(4.6+0.1) estimated by taking into account the factor 
2 uncertainty because of the inclusion of the NS atmosphere (Haensel 
2001)), and by considering the possible longer distance from the source 
to the earth d ~ 3.9fcpc (see, e.g. Pavlov et al 2004), respectively. The 
shadowed area stands for the possible M — R range of 1E1207.4-5209, 
which covers the strange matter and the normal neutron EOSs. 



m = 0.49(A/0.7)' (r^g/sJ 1 + Q.lh{A/Q.7Y Rl,,}j . (11) 

The quiescent spectrum of the transient type-I X-ray bursting 
NS Aql X-1 is measured with Chandra/ACIS-S, and the best-fit 
value of the apparent radius i?oo=13.4t4 (d/5fcpc) km is in- 
ferred by assuming a pure hydrogen atmosphere with the tem- 
peratures ranging from 145-168 eV, plus a power-law compo- 
nent. The distance from the observer to the source is 5 kpc as 
a fiducial value because the current uncertainties allow for a 
distance from 4 kpc to 6.5 kpc (Rutledge et al. 2001). There- 
fore, including the uncertainty on source distance the lower and 
upper limits on i?ooare 7.5 km ((13.4 - 4)* 4/5) and 23.9 km 
((13.4+5)*6.5/5), respectively. Moreover, the single kHz QPO 
frequency 1040 Hz has been detected in Aql X-1 by RXTE (e.g. 
van der Klis 2000, 2006), so that according to Eqs. (6) and Q, 
one can get the constraints A > 0.56 and m < 2.1. In Fig- 
urelD the ranges of M and R of NS in Aql X-1 are confined in 
the shadowed area where the meanings of the boundaries are 
indicated in the figure. We find that the shadowed area has a 
loose constraint on the NS EOSs, and it covers the represen- 
tative chosen EOSs of the strange matter, the normal neutrons 
and the pion condensate in the star core. Moreover, it is worth 
noticing that the A — 1.0 parabola divides the M — R diagram 
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Figure 2. Similar to Figure[T] The mass and radius of NS in Aql X-1 are 
constrained in tlie sliaded area. Rl (R2) represents Roo = 7.5(23.9) 
(km), wiiere the lower and upper hmits of iJooOn account of the uncer- 
tainty by the source distance are 7.5 km ((13.4 - 4)* 4/5) and 23.9 km 
((13.4+5)*6.5/5), respectively. The two concave-up parabolas stand for 
A = 0.56 and A = 1.0, respectively (as marked). The horizontal line 
m = 2.1 represents M=2.1M0. 



Figure 3. Similar to Figures □andU] The mass and radius of the NS in EXO 
0748-676 are constrained by the line z = 0.35 and two horizontal lines 
representing one (three) solar mass(es). Cottam et al. (2002) declaimed the 
gravitational redshift z = 0.35 with small error, but this samll error is not 
given quantitatively. As a possible comparison, we plot the hnes with the 
redshift 2 = 0.3 and z=0.4. 



into two parts (Figure. |2), and tiie strange star EOSs are rele- 
vant only when ^ > 1.0 

2.3 Constraining M and R with A and z: the case of EXO 
0748-676 

If the gravitational redshift z and the parameter A are known, both 
M and 7? can be calculated by the following two formulae, derived 
from Eqs.l2l( and (O, 

R(, ^ lA'if'\z){A/Q.7)-\ (12) 

m = lA'if^'^{z){A/Q.7)-^ . (13) 

In case that only a single kHz QPO frequency is detected, we can 
place a lower limit on the parameter A from Eq.([6l(. Following 
Eqs.(T2j and J13l >. M and R are constrained by the following in- 
equalities 

i?6 < (1850(Hz)/;/A')/"''(^) , (14) 

m < (1850(Hz)///k)/^''^(2) . (15) 

A gravitational redshift z — 0.35 was detected in EXO 
0748-676 by Cottam et al. (2002), who identify the most signif- 
icant features with the Fe XXVI and XXV n=2-3 and O VIII 
n=l-2 transitions, all with a redshift ot z = 0.35, identical 
within small uncertainties for the respective transitions, how- 
ever this small error has not been given quantitatively by the 
authors. In order to present the possible influence of the red- 
shift variation on the i\f—i? region, we plot z = 0.3 and z — 0.4 
lines in Figure|3l A single kHz QPO frequency 695 Hz is detected 
by RXTE for the same source (e.g., Roman & van der Klis 2000; 
van der Klis 2000). Thus A > 0.38 and m < 3.2 are inferred from 
Eq.l|6j and Eq.Q. As shown in Figure[3l the A > 0.38 condition 
gives a very loose constraint, which has little help to constrain the 



NS mass into a domain of less than 3 M© because the detected 
QPO frequency of this source is abnormally lower than the typ- 
ical values 1000 Hz (see, van der Klis 2000, 2006). If the future 
detections of this source present the higher frequency, then the 
promising constrain condition would be improved. Moreover, 
since M=0.97 ± 0.24Mq was measured as a lower mass for 2 A 
1822 - 371 (Jonker et al 2003), we set 1.0 Mq and 3.0 Mq as 
the conservative lower and upper mass limits for NSs, respectively. 
The M - R constraint for EXO 0748-676 is then limited to the 
segment of the straight dashed line z — 0.35 between the two hori- 
zontal lines m = 1.0 and jn = 3.0, where the EOSs for the strange 
matter, normal neutron and pion condensate in the star core are all 
permitted. 



3 SUMMARY AND DISCUSSION 

In this paper, we discuss the possibility of using any two of the 
three M ~ R relationships to constrain NS mass and radius. The 
three M ~ R relationships include the "apparent radius" (Roo), 
the gravitational redshift (z) and the averaged mass density (A^). 
We applied the method to three sample sources, 1E1207.4-5209, 
Aql X-1 and EXO 0748-676. In principle, if all the above three 
Al — R relations are well measured in one source, one can pose 
tight constraints on the ranges of M and R, so that the potential NS 
EOSs would be inferred. Unfortunately, this is not achieved so far. 
Furthermore, some M — R relations are incomplete or have large 
errors, so that additional astrophysical arguments are needed, and 
one can derive some inferred lower or upper limits of the M — R 
relations (e.g. Eqs.(6l and |7]). In binary systems, if the NS mass is 
known, any one of the three M — R relation measurements would 
be sufficient to derive its radius. 

For the compact star in 1E1207.4-5209, the estimates of its 
mass and radius come from the possible lower limit of the calcu- 
lated apparent radius of about 4.6 km (Bignami et al. 2004). An ex- 
traordinary low mass is inferred, which is even less than the mea- 
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sured lower limit of NS mass 0.97 ± O.24M0 in 2A 1822 - 371 
(Jonker et al. 2003). In such a case, the existence of exotic matter 
inside the star is hinted under the assumptions of blackbody spec- 
trum and 2kpc distance. 

The stellar radius can be also estimated by the motion-induced 
Doppler-broadening of absorption lines, when the emission region 
moves towards and away from the observer as the star rotates. 
According to this, Villarreal and Strohmayer (2004) have studied 
the source EXO 0748-676 with the measured gravitational redshift 
z — 0.35 (Cottam et al. 2002), and obtained a radius of about 1 1.5 
km and a stellar mass of about 1.75 M0 . A remark is that the emis- 
sion latitude and the orientation of the rotation axis both affect the 
result, so that the method is model-dependent. In addition, from the 
recent further study, Chang et al. (2005) show that the intrinsically 
broad line profile prohibits any meaningful constraint on the NS ra- 
dius if the 45 Hz burst oscillation seen by Villarreal & Strohmayer 
(2004) is the spin frequency. 

It is worth to mention the measured apparent radius of iso- 
lated NS in RX J1856.5-3754 because the distance of this source 
to earth is well determined to be about 117 pc (see, e.g. Walter 
& Lattimer 2002), with the perfect black-body spectrum of this 
source (Drake et al. 2002; Burwitz et al. 2003, 2001) the conserva- 
tive lower limit of apparent radius has been implied to be 16.5 km 
(d/1 17 pc) (Truemper et al. 2004). This corresponds to the "true" ra- 
dius of 14 km for a 1.4 M© neutron star, indicating a stiff equation 
of state at high densities, which excludes the quark stars or even 
the NSs with the quark cores (Truemper et al. 2004). Moreover, 
without knowledge of the spectrum line in RX J 1856.5-3754 un- 
til now (see e.g. Drake et al. 2002), we cannot constrain its M and 
R independently (Truemper et al. 2005). Once again, it is noted 
that the above conclusion comes from the black-body spectrum 
assumption, and the apparent radius may confront a modifica- 
tion by more than a factor of 2 if using the atmosphere model 
of photon spectrum (Haensel 2001). 

The next generation X-ray observatories with enhanced spec- 
tral and timing capabilities would greatly improve the "measure- 
ments" of the three M ~ R relations discussed in this paper. This 
would strengthen our confidence in evaluating the super-dense nu- 
clear compositions of compact stars. Theoretically, if both M and 
i? of a compact star are accurately measured, its EOS and nu- 
clear matter compositions would be explicitly revealed. This in 
turn sheds light on the unknown astrophysical processes during the 
supernova explosion of its progenitor star (e.g. Podsiadlowski et 
al. 2005). On the other hand, the regime of non-Newtonian strong 
gravitational fields would be also revealed with the determination 
of the NS mass and radius, where the direct tests of Einstein's the- 
ory of general relativity will be possible in the near future (e.g. van 
der Klis 2000, 2006). 
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